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PREFACE 


The  purpose  of  this  study  was  to  develop  a  method  for 
designing  straight-tube  evaporators  for  use  In  spacecraft  cooling 
systems.  Although  this  study  addressed  the  specific  problems 
associated  with  designing  an  evaporator  for  0-g  operation,  the 
method  of  design  is  applicable  to  all  straight-tube  evaporators. 

While  doing  this  study  and  writing  this  thesis  I  have  had  a 
great  deal  of  help  from  others.  I  would  like  to  especially  thank 
my  thesis  advisor,  Dr  J.  E.  Hitchcock  for  his  assistance  through 
all  aspects  of  this  study.  Furthermore,  I  would  like  to  thank 
Larry  Coulthart  of  the  Flight  Dynamics  Laboratory  for  his 
assistance  during  the  experimental  portion  of  this  study.  And 
finally,  I  would  like  to  thank  Julie  Roberts  for  her  typing  help 
and  patience  without  which  I  never  would  have  got  this  thesis  put 
together. 


ROBERT  E.  EASTMAN 
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ABSTRACT 


This  study  Involved  researching  the  design  of  straight-tube 
evaporators  for  spacecraft,  cooling  systems.  The  forced- 
convective-bolling  phenomenon  was  studied  to  determine  the 
characteristics  of  boiling  flow  and  how  gravitational  acceleration 
affects  this  process.  In  addition  to  this  theoretical  study,  an 
experimental  evaporator  was  built  to  investigate  the  operation  of 
straight-tube  evaporators.  It  was  found  that  pressure  fluctuations 
caused  by  slug  flow  may  hamper  efforts  to  build  two-phase,  cooling 
systems.  To  overcome  this  problem,  the  use  of  upstream  orifices 
and  twisted-tape  inserts  to  stabilize  the  flow  were  also 
investigated.  Also,  a  method  for  sizing  an  evaporator,  based  on 
the  best  available  method  of  calculating  the  local,  heat-transfer 
coefficients,  was  developed.  And  finally,  all  of  this  information 
was  used  to  develop  a  computer  program  to  design  straight-tube 
evaporators. 


I  INTRODUCTION 


In  the  past,  spacecraft  waste-heat  loads  have  been  relatively 
small.  For  this  reason  their  cooling  systems  have  been  rather 
simple  and,  in  many  cases,  fairly  inefficient.  Where  the  heat 
loads  of  the  past  and  present  spacecraft  have  been  at  most  25 
kilowatts,  many  future  applications  call  for  waste-heat  loads  as 
high  as  a  megawatt^.  Present  technology  cannot  be  scaled  up  to 
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meet  these  needs  without  severe,  weight  penalties  .  Therefore, 
to  meet  the  cooling  requirements  of  future  spacecraft,  new  methods 
of  transporting  waste-heat  loads  need  to  be  developed. 

The  easy  method  of  transporting  waste  heat  is  by  simple 
conduction  through  the  structural  members  of  the  spacecraft.  This 
has  proven  to  be  an  extremely  reliable  method  for  very  small, 
waste-heat  loads.  Unfortunately,  the  thermal  structure  quickly 
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gets  very  heavy  for  anything  but  a  small,  waste-heat  load1  \ 

When  the  waste-heat  load  grew  too  large  to  use  conduction 
through  the  structure,  thermal-design  engineers  began  to  use 
liquid-cooling  loops.  This  system  which  is  similar  to  the  cooling 
system  in  a  car,  has  worked  very  well.  In  fact,  the  space  shuttle 
uses  a  one-phase-f luid  loop  to  transport  its  waste-heat  load  to  the 
radiators.  This  system  could,  theoretically,  be  used  to  handle 
extremely  large,  waste-heat  loads.  Unfortunately,  as  these  systems 
get  large,  they  get  very  heavy.  Not  only  is  the  actual  cooling 
loop  heavy,  but  the  weight  penalty  associated  with  the  pump  drive 
also  gets  extremely  large.  In  fact,  in  many  cases,  the  equipment 
necessary  to  power  the  pump  is  much  heavier  than  the  cooling  loop 
Itself.  Futhermore,  the  temperature  drop  between  the  source  and 


the  radiator  is  quite  large.  Since  all  the  heat  must  be  radiated 
to  space,  the  radiator  also  gets  very  big  and  heavy. 

The  final  method  presently  available  for  transporting  waste 
heat  Is  the  heat  pipe.  A  heat  pipe  Is  nothing  more  than  a  tube, 
capped  at  both  ends,  that  contains  a  small  amount  of  liquid  and  a 
wick.  When  heat  Is  added  to  one  end  (the  evaporator)  the  liquid 
bolls.  Because  of  the  Increase  In  vapor  pressure,  the  vapor  Is 
transported  to  the  other  end  (the  condenser)  where  it  condenses  as 
heat  is  taken  out.  The  fluid  Is  then  wicked  back  to  the  evaporator 
to  complete  the  cycle.  The  heat  pipe  is  an  extremely  efficient  and 
reliable  method  of  transporting  heat.  The  problem  with  heat  pipes, 
however.  Is  that  they  can  only  transport  a  limited  amount  of  heat 
for  a  limited  distance^.  It  Is  possible  to  put  together  many 
heat-  pipe  bundles  In  series  to  overcome  this  problem,  but,  again 
weight  becomes  a  problem^. 

Due  to  the  shortcomings  of  the  previous  methods.  It  has  been 
proposed  to  use  two-phase-fluid  loops^.  Like  a  one-phase-fluid 
loop,  these  systems  would  involve  pumping  a  fluid  around  a  loop. 
Unlike  one-phase  systems,  two-phase  systems  would  use  the  latent 
heat  of  vaporization  rather  than  the  sensible  heat  of  the 
circulating  fluid.  Since  the  latent  heat  Is  much  greater  than  the 
sensible  heat,  this  system  can  use  a  much-smaller  flow  rate.  In 
fact,  preliminary  calculations  have  shown  that  the  pump  power 
required  will  be  approximately  two  orders  of  magnitude  less  than  a 
comparable  one-phase  system.  Furthermore,  the  entire  system  is 
essentially  Isothermal  which  will  mean  much  smaller  and  lighter 
radiators.  This  system,  however,  has  many  design  problems  that 


must  be  worked  out  before  it  can  be  used.  Therefore,  many  research 
programs  have  been  started  to  solve  the  problems  associated  with 
two-phase  systems. 

This  study  deals  with  the  design  of  the  evaporator.  More 
specifically,  the  objects  of  this  program  are  to  review  the 
force-convective  boiling  process  to  determine  the  effect  that 
gravity  will  have  on  it,  to  set  up  an  experimental  evaporator  to 
investigate  the  general  performance  and  problems  of  straigh-tube 
evaporators,  to  investigate  methods  of  solving  the  problems 
associated  with  them,  to  review  the  methods  available  for 
calculating  local  heat-transfer  coefficients  to  develop  a  best 
possible  approach,  to  develop  a  method  for  designing  straight-tube 
evaporators,  and  to  develop  a  computer  program  to  design 
straight-tube  evaporators. 


II  REVIEW  OF  THE  FORCED-CONVECTIVE  BOILING  PROCESS 


In  order  to  understand  the  relationship  between  the 
forced-convective  boiling  process  and  gravitational  acceleration,  a 
thorough  study  of  the  boiling  process  was  conducted.  The 
forced-convective  boiling  process  is  explained  and  how 
gravitational  acceleration  affects  the  heat-transfer  process  is 
examined  in  this  section. 

GENERAL  DISCUSSION 

The  forced-convective  boiling  process  can  be  broken  down  into 
five  regions:  one-phase,  forced  convection;  subcooled  boiling; 
saturated,  nucleate  boiling;  two-phase,  forced  convection;  and  the 
liquid-  deficient  region.  This  process  is  shown  in  Figure  1  and 
the  variation  of  the  heat-transfer  coefficient  associated  with  each 
region  is  shown  in  Figure  2.  The  one-phase,  forced-convection 
region  is  the  initial  condition  normally  encountered  in  an 
evaporator.  It  is  possible,  however,  that  this  region  may  be 
absent.  In  this  region  heat  is  conducted  into  and  is  swept  away  by 
the  fluid. 

As  the  temperature  of  the  wall  increases  along  the  length,  it 
reaches  a  point  where  bubble  nucleatlon  can  begin.  This  is  the 
start  of  the  subcooled-boiling  region.  As  the  wall  temperature 
increases  further,  the  bubbles  grow  larger  and  begin  to  depart  the 
surface.  Since  the  bulk  of  the  fluid  is  still  subcooled,  the 
bubbles  condense  upon  leaving  the  warmer  boundary  layer.  Although 
the  net  amount  of  vapor  generated  in  this  region  is  very  small,  the 
heat-transfer  coefficient  is  very  large. 
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Figure  2:  Variation  of  the  heat- transfer  coefficient  along  a  tube  during 
forced-convective  boiling. 


When  the  bulk  fluid  temperature  reaches  saturation  conditions, 
saturated-nucleate  boiling  occurs.  In  this  region  normal  nucleate 
boiling  occurs  with  a  corresponding  Increase  in  vapor  quality. 

Over  this  region  the  heat-transfer  coefficient  stays  essentially 
constant. 

As  more  vapor  is  produced,  the  liquid  layer  on  the  surface 
becomes  thinner.  Eventually,  this  liquid  layer  becomes  thin 
enough  to  prevent  the  liquid  immediately  adjacent  to  the  wall  from 
superheating  to  the  point  of  allowing  bubble  nucleation.  At  this 
point  the  mechanism  of  heat  transfer  switches  from  boiling  to 
liquid-surface  evaporation.  Since  bubble  nucleation  is  completely 
suppressed  and  the  heat  is  carried  away  by  the  vapor  at  the  liquid/ 
vapor  interface,  this  region  is  called  the  two-phase,  forced- 
convection  region.  This  region  has  the  highest  heat-transfer 
coefficients. 

The  final  region  occurs  when  the  vapor  quality  and  velocity 
increases  to  a  point  where  all  of  the  liquid  becomes  entrained  in 
the  vapor.  At  this  point  the  walls  dry  up  and  the  heat-transfer 
coefficient  decreases  dramatically.  This  "bum  out"  point  must  be 
avoided  at  all  costs  in  evaporator  design.  Heat  transfer  in  the 
liquid-deficient  region  will  not  be  considered  in  this  study. 

In  the  following  text  two  different  boundary  conditions  will  be 
used  to  further  illustrate  the  boiling  process  in  each  region.  The 
first  boundary  condition  will  be  a  constant  heat  flux  at  the 
evaporator  wall.  The  second  boundary  condition  will  be  a  constant 
evaporator  wall  temperature.  In  reality  a  combination  of  both  will 
be  seen;  but,  for  most  design  problems,  one  of  the  two  boundary 


conditions  can  be  used  with  reasonable  accuracy. 

SUBCOOLED  BOILING 

Subcooled  boiling  occurs  when  the  liquid  adjacent  to  the  wall 

is  superheated  to  a  point  where  bubble  growth  can  occur.  The 

heat-transfer  coefficient  associated  with  this  region  is  much 

larger  than  in  the  one-phase  region.  This  heat-transfer 

coefficient  increases  until  saturated-nucleate  boiling  is  achieved. 

The  heat-transfer  coefficient  begins  to  Increase  even  before 

bubble  nucleation  can  be  detected.  A  photographic  study  by  Jiji 
r  3 1 

and  Clark1  demonstrated  this  by  showing  that  a  considerably- 
higher  surface  heat  flux  is  required  to  permit  visual  detection  of 
bubbles  than  to  cause  an  improvement  in  the  heat-transfer 
coefficient.  This  increase  is  caused  by  hot  jets  of  fluid  from  the 
boundary  layer  being  projected  into  the  cooler  core  of  the  flow. 

For  the  case  of  a  constant  heat  flux,  the  onset  of  subcooled 
boiling  occurs  in  the  following  manner.  In  the  one-phase,  forced- 
convection  region  the  thermal-boundary  layer  next  to  the  surface 
heats  up  as  more  heat  is  absorbed.  Eventually  a  point  is  reached 
where  bubble  nucleation  can  occur.  A  bubble  will  grow  when  the 
pressure  inside  the  bubble  is  larger  than  the  pressure  outside. 

The  equilibrium  condition  for  the  spherical,  vapor  bubble  of  radius 


where  P is  the  pressure  of  the  vapor,  is  the  pressure  of  the 
liquid,  and  a  is  the  surface  tension  oi  the  fluid.  The  liquid 
superheat  required  to  maintain  this  equilibrium  can  be  related  to 
the  pressure  difference  with  the  use  of  the  Clausius-Clapeyron 
thermodynamic-equilibrium  relation: 


P 

v 


AT 


sat 


h.  AT  p  h. 

. fS  aat^v  fg 


T«.t 


(II-2) 


sat 


where  AT  _  is  the  superheat  of  the  liquid  at  the  wall  (T  -T  _). 
sat  w  sat 

On  any  surface  there  are  minor  imperfections  that  contain  a 
trace  of  vapor  or  noncondensable  gas.  These  imperfections  deter¬ 
mine  the  superheat  required  for  bubble  nucleation.  Therefore,  by 
combining  equations  (II-l)  and  (II-2): 
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it  can  be  seen  that  the  larger  the  size  of  the  Imperfection  the 
lower  the  superheat  that  is  required  for  nucleation. 
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Hsu  postualted  that  bubble  nuclei  on  a  heated  wall  will  grow 
only  If  the  lowest  temperature  on  the  bubble  surface  is  greater 
than  the  superheat  required  for  growth.  Earlier  work  by  this 
investigator^  confirmed  this.  For  a  bubble  of  radius  rc  to  grow 
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in  a  uniformly-heated  fluid  of  temperature  T^1 
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For  a  heated  surface  the  temperature  gradient  can  be  approximated 
in  a  linear  form:^ 


Tt(y>  -  T 


W 


(II-5) 


where  T£(y)  is  the  local,  liquid  temperature  and  y  is  the  distance 
away  from  the  wall.  For  a  nucleus  of  radius  r to  begin  to  grow 
equations  (II-4)  and  (II-5)  must  be  tangent^.  Therefore: 
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Bergles  and  Rohsenow^  obtained  a  graphical  solution  of  equations 
(II-4)  and  (II-5)  subject  to  equations  (II-6)  for  water  over  the 
pressure  range  of  15-2000  psia: 


W  -  15-6P1'156(4TBat) 


(2.3/P°- 02311 ) 


( II— 7) 


Davis  and  Anderson 


[7] 


developed  an  analytical  solution  of 


equations  (II-4),  (II-5),  and  (11-6): 


which  works  well  for  water. 


T81 

Frost  and  Dzakowic1  J  proceeded  in  the  same  manner  as  Davis  and 

Anderson;  but,  they  assumed  that  nucleation  occured  when  T^(y) 

2 

equaled  at  a  distance  of  r(Pr)  rather  than  rc>  Thus: 
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This  equation  compares  favorably  with  experimental  data  for  a 
variety  of  fluids. 

As  the  temperature  increases  further,  bubbles  form  on  the 

surface,  but  condense  before  they  detach.  During  this  phase  jets 

of  warmer  liquid  continue  to  flow  into  the  cooler  core  by  thermo- 

capillarity.  Thermocapillarity  is  caused  by  the  fact  that  the 

surface  tension  varies  around  a  bubble  growing  in  a  temperature 

gradient.  The  resultant  force  pulls  the  interfacial  liquid  film 

outward  (from  hot  to  cold)  which  then  drags  adjacent  warmer  liquid 

with  it.  Cooler  fluid  then  rushes  in  to  replace  the  warmer  fluid. 

[9] 

Brown  suggested  that  this  phenomenon  can  transfer  a  considerable 
fraction  of  the  heat  flux.  The  turbulence  induced  by  thermo¬ 
capillarity  explains  why  the  average  temperature  of  the  boundary 
layer  is  below  saturation.  Thermocapillarity  also  causes  the 
temperature  of  the  interface  between  the  one-phase  core  and  the 
bubble-laden,  boundary  layer  to  fall  with  increasing  heat  flux 
because  of  Increased  turbulence. 

In  the  final  stage  of  subcooled  boiling,  when  the  temperature 
approaches  saturation,  the  bubbles  grow  and  detach  from  the 


11 


surface.  However,  they  still  condense  as  they  migrate  into  the 
cooler  core  of  the  flow.  The  major  forces  affecting  the  flow 
involved  in  subcooled  boiling  are  the  force  exerted  by  the  fluid 
flowing  through  the  tube  and  the  thermocapillary  force  at  the 
thermal -boundary-layer  interface.  Neither  of  these  forces  are 
effected  by  gravitational  acceleration.  Furthermore,  gravitational 
acceleration  has  no  effect  on  bubble  nucleation.  For  these  reasons 
the  heat-transfer  coefficient  measured  in  a  1-g  environment  will  be 
the  same  as  that  measured  in  a  0-g  environment  for  subcooled, 
force-convective  boiling. 

In  actual  practice  the  equations  presented  in  this  section  may 
predict  lower  than  what  is  observed,  especially  at  lower 
pressures  and  low-heat  fluxes.  This  is  due  to  the  fact  that  the 
equations  assumed  a  wide  range  of  activation  sites.  On  metal  or 
glass  tubes  there  may  not  be  large  enough  nucleation  sites  for 
boiling  to  occur  even  though  the  equations  predict  it. 

For  the  case  where  there  is  a  constant  wall  temperature, 
subcooled  boiling  will  begin  immediately.  However,  the  boiling 
process,  discussed  above,  will  be  the  same. 

SATURATED-NUCLEATE  BOILING 

When  the  bulk-fluid  temperature  reaches  saturation  conditions 
the  vapor  bubbles  no  longer  condense.  Instead  they  continue  to 
grow  and  coalesce  forming  even  larger  bubbles.  Therefore,  the 
vapor  quality  begins  to  increase. 

In  this  region  the  heat-transfer  mechanism  is  very  similar  to 
nucleate-pool  boiling.  In  fact,  the  flow  rate  of  the  fluid  has 


very  little  influence  on  the  heat-transfer  coefficient^. 

Assuming  that  the  heat  transfer  is  caused  by  two  additive  effects: 
forced-convection  and  nucleate  boiling;  the  forced-convection 
component  can  be  ignored.  Therefore,  an  equation  for  the 
heat-transfer  coefficient  in  nucleate,  pool  boiling  can  be  used  to 
calculate  the  heat-transfer  coefficient  for  forced-convection, 
saturated,  nucleate  boiling.  This  has  been  done  successfully  by 
several  investigators^^’  However,  one  must  be  careful  in 

this  approach  because  subcooling  greatly  affects  pool  boiling, 
while  having  some  subcooling  (the  average  temperature  is  the 
saturation  temperature)  in  the  liquid  core  for  forced  convection 
has  no  effect. 

The  reasons  for  the  large  increase  in  the  heat-transfer  coeffi¬ 
cients  for  nucleate  boiling  have  been  debated  for  years  with  no 
concise  agreement^.  The  main  mechanism,  however,  seems  to  be 
the  turbulence  caused  by  bubble  generation  and  departure.  Before 
the  bubble  grows,  a  thermal-boundary  layer  is  built  up.  When  the 
temperature  near  the  surface  gets  high  enough  for  bubble  growth, 
the  bubble  grows  outward,  pushing  the  boundary  layer  with  it. 
Furthermore,  when  the  bubble  detaches  from  the  surface  it  then 
carries  this  thermal-boundary  layer  into  the  bulk  of  the  fluid, 
thus,  directly  transporting  energy  outward.  Cooler  fluid  then 
rushes  back  to  the  surface  to  complete  the  cycle. 

The  major  factors  affecting  heat  transfer  during  saturated- 
nucleate  boiling^ ^  are  the  heat  flux  and  the  pressure.  The  flow 
rate,  quality,  and  core  subcooling  have  little  effect  in  this 
region.  Furthermore  the  heat-transfer  coefficient  remains  constant 
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through  out  the  region  because  of  the  constant,  bulk  temperature 
(equal  to  the  saturation  temperature).  And,  the  boundary 
conditions  at  the  wall  have  no  effect. 

In  this  region  the  major  forces  acting  on  the  fluid^  are 
surface  tension,  Inertia,  bouyancy,  drag,  and  the  Bernoulli  effect. 
Surface  tension  holds  the  bubble  to  the  wall  while  the  Inertia  of 
the  liquid  surrounding  a  growing  bubble  attempts  to  pull  it  from 
the  wall.  The  drag  force  caused  by  the  fluid  flowing  around  the 
bubble  attempts  to  shear  It  from  the  wall.  The  bouyancy  force 
pulls  the  bubble  upward,  pulling  It  from  the  wall  or  pushes  It  onto 
the  wall  depending  on  the  surface  orientation.  Upon  departure  from 
the  surface,  bouyancy  makes  the  bubble  rise  while  the  Bernoulli 
effect  created  by  the  velocity  gradient  within  the  tube  causes  the 
bubbles  to  migrate  to  the  center.  As  long  as  the  flow  Is  not 
stratified,  the  differences  between  the  heat-transfer  coefficient 
In  1-g  and  0-g  for  a  round  tube  should  be  small.  This  Is  because, 
in  1-g,  gravity  helps  on  the  bottom  of  the  tube  by  helping  to 
remove  the  vapor  but  hinders  on  the  top  by  keeping  the  vapor  near 
the  surface.  The  average  effect  should,  therefore,  be  similar  to 
0-g  conditions. 

TWO-PHASE,  FORCED  CONVECTION 

As  the  vapor  quality  Increases  the  vapor  migrates  towards  the 
center  creating,  at  first,  a  slug-flow  pattern  which  Is  followed  by 
an  annular-flow  pattern.  A  slug-flow  pattern  Is  where  large 
bubbles  or  "slugs"  of  vapor  are  separated  from  one  another  by  thin. 
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liquid  regions.  As  the  vapor  quality  increases  the  liquid  regions, 

separating  the  slugs,  disappear  leaving  an  annular-flow  pattern 

where  the  vapor  core  is  separated  from  the  wall  by  a  thin,  liquid 

layer.  Flow  patterns  and  the  conditions  for  their  appearance  is 

discussed  extensively  in  Reference^ 

Eventually,  a  point  is  reached  where  the  liquid  layer  adjacent 

to  the  wall  is  thin  enough  that  all  of  the  heat  flux  is  conducted 

through  the  liquid  without  nucleate  boiling.  The  mechanism  of  heat 

transfer  then  switches  from  nucleate  boiling  to  evaporation  at  the 

liquid/vapor  interface.  This  regime,  therefore,  is  called  the 

two-phase,  forced-convection  region. 

In  order  to  suppress  bubble  nucleation,  the  effective  thermal 

conductivity  of  the  liquid  layer  must  be  high  enough  so  that  the 

liquid  cannot  be  superheated  to  a  point  where  bubble  growth  occurs. 

The  equations  used  earlier  to  predict  the  onset  of  subcooled, 

nucleate  boiling  are  applicable  here  too.  As  mentioned  earlier, 
f  81 

Frost  and  Dzakowic1  1  predicted  the  bubble  growth  would  begin  when: 

4ora  •  (~pH 
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Combining  this  equation  with  Newton's  law  of  cooling: 
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(11-10) 


where  ^Q^ppQ  is  the  heat-transfer  coefficient  for  the  onset  of 
two-phase,  forced  convection,  and  q^g  Is  the  heat  flux  for  the 
suppression  of  nucleate  boiling.  Setting  <1qNq  * 
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for  for  the  onset  of  two-phase,  forced  convec¬ 

tion  can  be  found. 
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This  equation  for  h0TpFC  is  the  point  where  bubble  nucleation 
is  suppressed.  As  in  subcooled  boiling,  this  number  may  be  lower 
than  what  is  actually  found  in  practice  due  to  the  absence  of  large 
nucleation  sites. 

In  this  two-phase,  forced-convection  region  the  heat  is  trans¬ 
ferred  across  the  liquid  layer  and  evaporates  liquid  at  the  liquid/ 
vapor  interface.  Therefore,  the  heat-transfer  coefficient  is  very 
dependent  on  the  thickness  of  the  liquid  layer.  This  thickness  is 
governed  by  the  quality.  Therefore,  h^pp^  is  very  dependent  on 
quality.  Furthermore,  as  the  thickness  of  the  liquid  layer 
decreases,  the  temperature  drop  across  the  liquid  layer  decreases. 
Because  of  this  extremely-hlgh,  heat-transfer  coefficients  can  be 
achieved. 

Another  factor  that  affects  the  heat-transfer  coefficient  is 
the  mass  flux^.  Higher  mass  fluxes  Increase  the  vapor  velocity 
which  increases  the  evaporation  at  the  liquid/vapor  Interface. 
Unlike  the  saturated-nucleate-boiling  regime,  the  heat  flux  has  no 
effect  on  the  heat-transfer  coefficient. 

Since  the  heat-transfer  coefficient  increases  along  the  tube 


with  increasing  quality,  the  wall  temperature  of  the  tube  will 
decrease  along  the  tube  for  the  case  of  a  constant  heat  flux.  For 
the  case  of  a  constant  wall  temperature,  the  heat  flux  will 


increase  along  the  tube. 


In  this  regime,  the  only  contribution  of  gravity  to  the  heat 
transfer  process  is  in  the  distribution  of  film  thickness  around 
the  circumference  of  the  tube.  However,  the  average  value  will  be 
the  same  as  the  constant  value  found  under  0-g  conditions^. 
Therefore,  heat  transfer  in  this  regime  should  be  relatively 
insensitive  to  gravity.  Furthermore,  all  analytical  treatments  of 
this  regime  have  used  a  constant,  fluid  thickness  around  the  peri¬ 
meter  (0-g  condition)  successfully  in  past  analyses^. 

LIQUID  DEFICIENT  REGION 

As  the  quality  increases  in  the  two-phase,  forced-convection 
region  the  liquid  layer  gets  thinner,  while  the  vapor  velocity 
increases.  Eventually  a  point  is  reached  where  the  liquid  becomes 
totally  entrained  in  the  vapor  core,  thus  causing  the  walls  to 
dryout.  When  dryout  conditions  are  reached,  the  heat-transfer 
coefficient  decreases  drastically.  With  a  constant  heat  flux,  this 
causes  the  wall  temperature  to  increase  tremendously.  Often  it 
will  increase  above  the  melting  point  of  the  evaporator  walls,  thus 
"burning  out"  the  evaporator.  For  the  case  of  a  constant  wall 
temperature,  the  heat  flux  decreases  drastically.  This  causes  the 
liquid  deficient  region  of  the  evaporator  to  be  essentially 
useless.  This  condition,  therefore,  must  be  avoided. 


Ill  FLOW-VISUALIZATION  EXPERIMENT 

In  order  to  observe  the  general  operation  of  and  problems 
associated  with  straight-tube  evaporators,  a  flow-visualization 
experiment  was  set  up.  The  experimental  apparatus  and  the  results 
of  the  experiment  are  discussed  below. 

EXPERIMENTAL  APPARATUS 

The  flow-visualization  loop  consisted  of:  an  evaporator,  a 
condenser/accumulator,  a  preheater,  a  pump,  a  flow  meter,  and  an 
adjustable  orifice.  This  set  up  is  shown  in  Figure  3.  The 
evaporator  consisted  of  a  long,  small -diameter,  glass  tube  and  will 
be  discussed  in  greater  detail  later.  The  condenser/accumulator 
consisted  of  a  copper  coil  (condenser)  submerged  in  a  bucket  of 
water  (accumlator).  The  steam  and  water  mixture  from  the  evapora¬ 
tor  flowed  through  the  copper  coil  and  emptied  Into  the  bottom  of 
the  bucket.  The  copper  coil  thus  served  as  a  condenser  and  as  a 
heater  for  the  water  supply  in  the  accumulator.  The  accumlator  was 
open  to  the  atmosphere.  The  preheater  was  a  hot  plate  beneath  the 
accumulator.  It  was  used  to  help  raise  the  temperature  of  the 
water  to  near  saturation.  The  flow  meter,  a  rotameter,  consisted 
of  a  small,  vertlcle,  glass  tube.  A  steel  and  a  rubber  ball  rose 
within  the  tube  with  an  Increase  in  the  flow  rate.  The  flow  meter 
also  contained  a  valve  that  could  be  used  to  regulate  the  flow. 

The  purpose  of  the  adjustable  orifice  will  be  discussed  later. 

Instrumentation 


Total  instrumentation  consisted  of  three  pressure  meters. 


Figure  3:  Experimental  apparatus 


three  thermocouples,  an  ammeter,  and  a  voltmeter.  Two  of  the 
pressure  meters  vere  set  up  to  measure  the  pressure  drop  across  the 
evaporator  and  across  the  adjustable  orifice.  The  third  pressure 
gage  measured  the  pressure  of  the  flow  prior  to  entering  the 
orifice.  The  three  thermocouples  vere  located  before  and  after  the 
evaporator  and  in  the  accumulator.  The  ammeter  and  voltmeter  were 
used  to  measure  the  power  input  to  the  evaporator. 


Evaporator 

The  evaporator  consisted  of  a  long,  thin,  glass  tube.  Heat  was 
supplied  to  the  tube  by  nichrome  heating  wire  which  was  coiled 
around  the  tube.  The  heat  from  the  wire  was  conducted  through  the 
glass  tube  which  caused  the  water  to  boll  as  it  passed  through  the 
tube. 

The  purpose  of  experiment  was  to  visualize  the  boiling  process 
in  a  spacecraft  evaporator.  To  do  this  it  was  necessary  to 
minimize  the  influence  of  gravity.  On  earth,  gravity  tends  to 
dominate  two-phase  flow  through  tubes,  while  in  space,  surface 
tension  tends  to  dominate.  However,  in  very  small-diameter  tubes, 
surface  tension  dominates  over  gravity.  For  this  experiment,  this 
point  was  assumed  to  be  the  point  where  the  liquid  is  able  to  wick 
itself  completely  around  the  tube  so  that  a  stratified-flow  pattern 
cannot  exist.  From  earlier  work^  this  point  is: 
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Using  water  properties: 


f  2(5. 
Tm[(955 
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(955kg/m^)(9.8a/s2)J 


r  -  0.0035m  «  O.lUft  -  l/8in 

Therefore,  tubes  with  %  inch  diameters  or  less  were  used  for  the 
evaporator . 


RESULTS  OF  FLOW  VISUALIZATION 


The  flow  entered  the  evaporator  subcooled.  Almost  immediately 
very  small  bubbles  would  appear  at  the  evaporator  wall  (subcooled 
boiling).  These  bubbles  initially  didn't  grow;  but,  further  along 
the  tube  as  the  bulk  temperature  increased,  they  began  to  grow, 
depart  from  the  surface,  and  recondense  in  the  flow.  As  the  bulk 
temperature  of  the  flow  approached  saturation,  the  bubbles  began  to 
grow  larger  and  coalesce.  When  the  bulk  temperature  reached 
saturation  the  flow  pattern  changed  from  bubbly  flow  to  slug  flow. 
In  this  slug-flow  regime  the  liquid  layer  between  the  vapor  core 
and  the  wall  was  thin  enough  to  suppress  nucleate  boiling.  Still 
further  down  the  tube,  the  slugs  expanded  into  an  annular-flow 
pattern.  This  pattern  established  in  the  evaporator  was  always  the 
same  regardless  of  power  input  or  flow  rate.  The  only  difference 
being  whether  enough  power  was  being  supplied  to  reach  the  flow 
quality  necessary  for  annular  flow. 

Although  the  flow  pattern  established  in  the  evaporator  was  consist 
ent,  it  was  not  stable.  This  instability  was  caused  by  the  extreme 
pressure  oscillations  associated  with  slug  flow.  In  a  spacecraft 
evaporator  these  pressure  oscillations  cannot  be  allowed.  They 


would  cause  vibration,  disruption  of  flow  uniformity  In  parallel- 
flow  networks,  wear  on  pumps,  and  instabilities  In  the  condenser. 
Therefore,  before  straight- tube  evaporators  can  be  used,  a  means  of 
stabilizing  the  flow  must  be  devised. 

In  an  effort  to  solve  this  problem,  an  orifice  was  installed 
Immediately  In  front  of  the  evaporator  tube.  If  the  orifice  could 
be  sized  so  that  the  pressure  drop  through  the  evaporator  Is  much 
smaller  than  the  pressure  drop  through  the  orifice,  the  flow  should 
be  stable  In  front  of  the  evaporator.  This  would  eliminate  the 
problems  of  parallel-flow  Instability  and  pump  wear.  To  test  this 
theory,  an  adjustable  orifice  was  placed  ahead  of  the  evaporator. 
The  pressure  was  measured  at  point  1  (see  Figure  3),  and  the 
pressure  drop  was  measured  between  points  1  and  2  and  between  2  and 
3.  The  results  of  this  test  were  very  conclusive  -  this  method 
will  not  work!  The  orifice  In  this  experiment  was  varied  until  the 
pressure  drop  across  the  orifice  was  approximately  an  order  of 
magnitude  greater  than  the  pressure  drop  through  the  evaporator. 

And  still,  the  flow  oscillations  could  not  be  damped  out  prior  to 
the  orifice.  The  pressure  readings  are  shown  in  Table  I  for 
various  mass  flow  rates.  An  investigation  of  the  accumulator 
showed  the  flow  to  be  steady  to  the  pump,  therefore  the  only  cause 
of  these  instabilities  could  be  the  slug  flow.  Any  further 
increase  in  the  pressure  drop  across  the  orifice  was  deemed 
unacceptable  because  the  corresponding  increase  in  pump  power 
required  would  eliminate  the  major  advantage  of  the  two-phase 
systems  having  a  low  pump-power  requirement. 
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IV  TWO-PHASE,  HEAT-TRANSFER  EQUATIONS 


In  order  to  determine  a  method  of  calculating  the  heat-transfer 
coefficients  for  forced-convective  boiling,  a  thorough  review  of 
past  work  Is  done.  The  Information  obtained  In  this  section  is 
used  to  develop  a  better  method  In  Section  V. 

Many  empirical  equations  have  been  developed  to  predict  the 
heat-transfer  coefficients  in  forced-convectior-boiling  heat 
transfer.  Most  of  the  equations  were  developed  for  only  one  region 
of  boiling,  while  others  use  weight  factors  with  several  equations 
to  make  a  general  equation.  The  weight  given  to  a  certain  equation 
depends  on  the  local  flow  conditions.  This  section  describes  some 
of  the  newer  and/or  more  popular  equations. 

One  of  the  easiest  equations  to  use  is  that  of  McAdams  et 
alJ12^  They  developed  an  empirical  equation  for  the  heat  flux  for 
water  during  subcooled,  nucleate  boiling: 

q  .  c'AT  3-86  ( IV-1) 

where  C'  is  a  dimensionless  constant  that  ranges  from  0.190  to 
0.074  depending  on  the  extent  of  degassing  of  the  water  and  T$flt 
is  the  difference  between  the  wall  temperature  and  the  saturation 
temperature.  They  found  q  to  be  independent  of  water  velocity, 
pressure,  degree  of  subcooling,  and  diameter.  The  equation  is  good 
for  pressures  of  30  psia  to  90  psia.  The  main  problem  with  this 
equation  Is  that  It  was  developed  solely  for  water  and  is  not 
applicable  to  other  fluids. 

As  mentioned  previously,  pool-boiling  equations  have  been  used 


successfully  for  cslculating  Che  heac-transfer  coefficients  in  the 


subcooled  and  saturated-nucleate-boiling  regions.  Specif ically, 

Lavin  and  Young ^ ^  ^ ,  and  Myers  and  Katz^**^  have  successfully  used 
[131 

McNelly's  equation1  1  for  the  nucleate,  pool-boiling  heat-transfer 
coefficient  to  correlate  their  saturated,  forced-convection 
nucleate-boiling  data.  McNelly's  equation  is: 


.0.69 


(—) '  ( 

’  Cp2.V 

\  (h.  A 

.  / 

1  U'.v 

0.31,^0.31 


where  C  was  given  as  .225  by  McNelly  but  modified  to  .25  for  the 
inside  of  a  plain  tube  by  Lavin  and  Young.  And: 
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McNelly  developed  this  equation  by  first  looking  very  closely 
at  the  boiling  process,  and  then  doing  a  dimensional  analysis.  He 
also  used  the  heat  flux  as  the  driving  mechanism  rather  than  the 
temperature  difference.  It  should  be  noted  that  like  McAdams 
equation  the  diameter  has  no  effect  (divides  out)  and  the  mass  flux 
has  no  effect  (developed  for  pool  boiling) .  This  equation  is  good 
for  all  fluids. 

Lavin  and  Young  also  went  on  to  develop  an  equation  for  the 
two-phase,  forced-convection  region.  To  correlate  their  data  they 
develop  the  equation: 


where  X  is  the  Lockhart-Martinelli  parameter  for  turbulent 
liquid,  turbulent  vapor: 


(IV-8) 


Dengler  and  Addoms 
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found  that: 


(IV-9) 


where  h^0  is  the  heat-transfer  coefficient  for  liquid  alone  found 
by  setting  x-0.0  and  using  equation  (IV-5) .  Guerrieri  and 
Talty^5^  found  that: 


(IV-10) 


where  h^p  is  the  one-phase,  heat-transfer  coefficient.  Bennett  et 

aifl6]  jcouncj  that: 


(IV-11) 


All  of  the  above  correlations  were  compared  by  Chen1  .  He 
found  that  Dengler  and  Addom's  equation  predicted  values  too  large 
when  compared  with  other  experimenters  data.  Guerrieri  and  Talty's 
equation  predicted  lower  coefficients  than  Dengler  and  Addom's  but 
still  lacked  accuracy.  And  Bennet's  predicts  h^p  well  for  water 
but  not  for  other  fluids.  For  these  reasons  Chen  developed  his  own 
method  which  covers  both  the  saturated,  nucleate-boiling  region  and 
the  two-phase,  forced-convection  region.  He  assumed  that  the  two 
mechanisms,  nucleate  boiling  and  forced  convection,  are  present  to 
some  extent  in  both  regions.  Therefore,  he  devised  a  correlation 
using  superposition: 


hTP  *  hNB  +  hFC 


(IV-12) 


Where  h^  is  the  contribution  due  to  nucleate  boiling  and  hpC  is 
the  contribution  due  to  forced  convection.  A  Dittus-Boelter  type 
equation  was  used  for  hpg. 
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Chen  argued  that  kpp  and  Pr^p  can  be  equated  with  the  liquid 
values.  Furthermore,  he  defined  a  parameter  F  as: 
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The  only  unknown  in  then  is  F.  However,  since  F  is  a  flow 
parameter  it  can  be  related  to  Chen  used  Foster  and  Zubers' 

analysis  of  nucleate  boiling  for  h^^a 
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He  also  defined  a  nucleate-boiling,  suppression  factor,  S,  which 
relates  the  mean  superheat,  AT,  to  the  wall  superheat,  At 


using  the  Clausius-Clapeyron  equation: 
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(IV- 17) 


Combining  equations:  (IV-15),  (IV-16)  and  (IV-17);  h^  equals: 
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ATS  (IV-18) 


The  Parameters  S  and  F  are  graphed  in  Figures  4  and  5.  The  Chen 

equation  has  been  correlated  with  the  data  of  many  investigators 

with  85Z  of  the  data  being  within  -  15. IX.  The  problem  with  the 

Chen  equation  is  that  it  is  not  easy  to  program  on  a  computer  due 

to  the  use  of  charts  to  find  S  and  F. 

A  final  method  for  calculating  h^p  was  proposed  recently  by 
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Bjorge,  Hall,  and  Rohsenow.  1  They  also  used  a  superposition 
formula  where: 


^TP  *  qFC  +  qNB  qIB 


(IV-20) 


for  the  high  quality  region  (x  >  0.05),  and 


(IV-21) 


for  the  low  quality  region,  q^  is  the  heat  flux  at  the  incipient 

3 

boiling  point.  Since  q^g  and  qjg  are  proportioned  to  T  : 


*  ^FC  + 


(IV-22) 


for  the  high  quality  region,  and  for  the  low  quality  region: 


Figure  S:  Suppression  factor  S  versus  the  Reynolds  number 


For  the  high  quality  region,  the  equation  developed  by  Traviss  et 
al^*^.  and  modified  by  Hall^^  was  used  for  q^: 


*FC 


<V!i i)kt 
V 


AT 


sat 


(IV-24) 


where : 


F(Xm) 


Kit  * 


(IV-25) 


O-  .  Gd(l-x)  (IV-6) 

Rei  v4 

F2  "  5Pr£+5 In ( 1+5  Pr£ ) +2 . 5 In ( 0 . 0  0  3 lRe ^ • 812 )  (Re>1125)  <lV-26a) 
F2  -  5Pr£+51n[l+Pr£(0.096i»Re4°*585-l)]  (50<Re<1125)  (IV-26b) 


F-  ■  O.OTOTPr.Re 


0.5 


(Re >50 )  <IV-26c> 


The  equation  vas  developed  by  Travlss  et  al.  as  a  correlation  for 
condensation  in  annular  flov.  Since  tvo-phase ,  forced-convection 
evaporation  is  the  opposite  of  forced-convection  condensation,  this 
equation  should  work  for  tvo-phase  evaporation.  Hall  found  this  to 
be  true,  with  a  slight  modification  of  the  coefficient. 

For  the  subcooled  and  lov-quality  region,  Colburn's  equation 
was  used: 


(IV-27) 


where  the  subscript  f  means  that  the  property  should  be  evaluated 
at  the  film  temperature,  (Tw  +  T^)/2,  the  subscript  b  refers  to  the 
bulk  temperature,  and  qpc  can  be  found  with  equation  (IV-3) 

The  nucleate  boiling  heat  flux  for  both  regions  is  calculated 
by  using  the  Mikic-Rohsenov  correlation: 


B 


1/2  17/8.  19/8  1/8 
p£  pv 


AT 


sat 


(IV-28) 
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where  is  1.89X10  in  SI  units  for  forced-convection  boiling  of 
water. 

For  the  wall  superheat  at  the  incipience  of  boiling,  they  used 


the  equation  derived  by  Davis  and  Anderson  for  the  high  quality 
region: 


In  the  subcooled  and  low  quality  region,  this  equation  Is  modified 
to  account  for  the  fact  that  there  may  not  be  large  nucleatlon 
sites: 


ATsat,IB  *  l-N^rs  “  RATsc)  ^rtang>rmax^  (IV“30) 

where  rtang  1*  the  radius  of  tangency  between  equation  II-4  and 
II-5  in  the  subcooled-boiling  section  or: 

UgTsatvfg 

rtang  “  bfgATgat4lB  (IV-31) 

Furthermore: 


r 


faTTJl  h. 


sat  fg  FC 


(IV-32) 


and 


N 


k„ 


(IV-33) 


Where  r  *  10 
max 

It  should  be  noted  that,  in  the  method  devised  by  Bjorge,  Hall 


and  Bohsenow,  hyC  for  the  high-quality  region  is  completely 
different  from  h__  in  the  lov-quallty  and  subcooled  regions.  In 


the  low-quality  and  aubcooled  region  represents  the  one-phase, 
forced-convection,  heat-transfer  coefficient,  equation  (IV-27). 
While  in  the  high-quality  region,  h^  represents  the  heat  transfer 
coefficient  for  forced  convection  with  a  phase  change  at  the 
boundary,  equation  (IV-24) . 


V  PROPOSED  METHOD  FOR  DESIGN  OF  STRAIGHT-TUBE  EVAPORATORS 


In  this  section  the  Information  that  was  presented  In  the 
preceedlng  three  sections  Is  combined  to  develop  a  new  approach  to 
designing  straight-tube  evaporators  and  a  new  method  for  sizing  the 
length  of  the  evaporator  tube  Is  developed.  Next  a  new  way  of 
calculating  the  local  heat-transfer  coefficients  and  how  to 
stabilize  the  flow  is  presented.  This  is  followed  by  a  way  to 
calculate  the  pressure  drop.  And  finally  all  of  this  information 
was  used  to  develop  a  computer  design  program  which  Is  also 
explained  in  this  section. 

DESIGN  METHOD 

Presently,  the  length  of  the  tube  necessary  for  a  certain  size 
of  evaporator  Is  calculated  by  treating  the  one-phase  and  two-phase 
sections  separately.  First  the  length  of  the  one-phase  section  is 
calculated  by  using  the  one-phase  heat-transfer  coefficient,  the 
energy  It  takes  to  reach  saturation  conditions,  and  the  desired 
wall  temperature.  This,  of  course,  gives  a  certain  heat  flux  in 
the  one-phase  section. 

After  the  calculations  are  finished  for  the  one-phase  section, 
the  length  of  the  two-phase  section  is  calculated  in  the  same  way, 
using  appropriate  two-phase  equations.  For  this  calculation  it  is 
assumed  that  boiling  begins  when  the  bulk  temperature  reaches 
saturation.  The  total  length  can  then  be  calculated  by  adding  the 
lengths  and  multiplying  by  a  safety  factor. 

The  problem  with  this  approach  is  that  by  breaking  up  the  tube 
into  two  separate  parts  the  heat  flux  is  considerably  different  in 


the  two  sections.  In  fact,  since  Q/A  is  proportional  to  h  and  h^p 
is  less  than  hTp,  Q/App  be  muc^  less  than  Q/A,^.  This, 

however,  will  not  be  the  case  in  an  actual  evaporator.  Instead, 
there  will  nearly  be  a  constant  heat  flux,  such  as  with  electronics 
cooling,  or  a  constant  wall  temperature,  such  as  with  exchanging 
heat  with  a  heat  pipe.  Furthermore,  with  a  constant,  wall 
temperature,  subcooled  boiling  will  occur  rather  than  one-phase 
heat  transfer.  With  a  constant  heat  flux  the  wall  temperature  will 
vary  with  the  heat-transfer  coefficient,  h,  while  with  a  constant 
wall  temperature  the  local  heat  flux  will  vary  with  h. 

Because  of  the  inaccuracies  of  the  past  approach,  the  following 
two  methods  of  "sizing"  an  evaporator  have  been  developed.  These 
methods  are  depicted  in  Figures  6  and  7  for  a  constant  heat  flux, 
and  for  a  constant  wall  temperature,  respectively.  For  a  constant 
heat  flux,  first  calculate  the  incipient-boiling  point.  This  is 
the  point  where  the  wall  temperature  is  high  enough  above  the 


saturation  temperature  that  boiling  can  begin.  Next,  calculate  h 


OP 


and  hTp»  h^,p  will  vary  along  the  length  of  the  two-phase  section. 


Using  the  allowable  At  ■  T  -T  ,  calculate  the  length  using 

S3l  W  S3t 


Qtotal’  avera8e  ^XP’  ^  an<*  ^sat'  *s  t*ie  length  of  tube 


necessary  for  the  evaporator.  The  length  of  the  subcooled  section 
is  the  length  necessary  to  get  ATtQta^-  ATgc  ■  ATJB,  where  ATtQta^ 


Tv  -  V  *T.c  •  T8at  -  V  *nd  *TIB  '  TIB  -  Tsat*  Fro"  the 


incipient  boiling  point  to  the  saturation  point,  subcooled  boiling 


will  take  place  and  hTp  must  be  used.  Prior  to  that  point  hgp  can 


be  used  to  determine  the  wall  temperature. 


For  a  constant  wall  temperature,  the  heat  flux  will  vary  along 
the  tube.  Since  the  wall  temperature  will  be  higher  than  the 
Incipient-boiling  point,  boiling  will  start  Immediately.  There¬ 
fore,  using  appropriate  two-phase  equations,  hTp  and  the  heat  flux 
can  be  calculated.  And  finally,  the  length  can  be  found. 


CALCULATION  OF  THE  HEAT-TRANSFER  COEFFICIENT; 

To  calculate  the  heat-transfer  coefficient,  the  superposition 
method  of  Rohsenow^®^  was  employed.  This  method  was  used  because 
of  its  success  in  the  past  and  its  adaptability  to  a  computer 
solution.  This  method  assumes  that  the  forced-convection  boiling, 
heat  flux  can  be  found  by  summing  the  heat  fluxes  for  pool  boiling 
and  forced-convection.  More  specifically,  for  subcooled  boiling 
flow  (x=0): 


I 

l 

(IV-23) 


and  for  saturated-boiling  flow  (x  >  0): 


^Tp 


fe2)] 


(IV-22) 


and  for  high-quality  flows  where  nucleate  boiling  is  suppressed 


FLUID  r 

CALCULATE: 

AT  incipient  boiling 

PROPERTIES 

' 

CALCULATE: 

h  one-ohase.  h  two-ohas  \ 

n  m  n  M  ■ 

CALCULATE: 

TOTAL  LENGTH 

h  two-phase 

CALCULATE: 

ONE  PHASE  LENGTH, 
LOCAL  TEMPERATURES 

_6:  Flow  chart  for  calculating  the  length  of  an  evaporator 
with  a  constant,  heat  flux. 


Figure  7:  Flow  chart  for  calculating  the  length  of  an  evaporator 
with  a  constant,  wall  temperature. 


Where  qPQ  for  the  subcooled  regions  is  the  one-phase, 
forced-convection  heat  flux,  equations  (IV-27)  and  (IV-2),  and  qpc 
is  the  forced-convection,  evaporative  heat  flux  equations  (IV-24), 
for  the  saturated  regions.  The  force-convective  heat  flux  for  the 
subcooled  region  can  be  calculated  by  Colburn's  equation: 


V  «  0.023' 


0.333 


(IV-27) 


wi  th : 


9FC  *  hFC(ATSat  +  ATsc} 


(V-2) 


qFC  for  the  high-quality  region  can  be  found  by  Traviss's 
equation: 


"■t0‘9prirl*ttV.- 

qFC  *  F„d  ATsat 


(IV-24) 


F(Xtt)  -  0.15 
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Gd(l-x) 


(IV-25) 
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F2  *  5PrJl+51n(l+5Prji)42.5ln(0.0031Re£0'®12)  (Re>1125)  (IV-26a) 
F2  *  5Pr£451n^l+Pr£(0.095J4Re£°-585_i)J  (50<Re<1125)  (IV-26b) 


F2  -  O.OTOTPr^  *' 
q^g  for  both  equations  can  be  found  by  using  McNelly's  equation: 


(JV-2) 


where  qTp  is  the  total  heat  flux  and: 

%B*  hBATBat 

And  finally,  AT$at  jg  can  be  found  with  Frost  and  Dzakowic’s 
equation: 


4T  .  ?fWwipr2 

sat, IB  k.h .  rr 

£■  fg 


(V-4) 


where  hpc  is  found  from  equation  (IV-27), 


COMPARISON  WITH  DATA 

The  method  outlined  above  was  compared  with  data  from  Cheng, 
et.al  J21^;  Latsh,  et.al  J22-^  and  Wright^23-^.  Cheng's  data  was  for 
subcooled  boiling  water  and  is  shown  In  Figure  8.  In  this  figure 
both  the  method  outlined  above  and  the  original  method  by  Bjorge, 
Hall  and  Rohsenow  is  shown.  Notice  that  there  is  very  little 


difference  between  the  two.  This  is  due  to  the  fact  that  both 
methods  are  the  same  except  for  the  equation  used  to  calculate 
qNB.  For  water  both  these  equations  give  essentially  the  same 
answers.  However,  the  equation  used  by  Bjorge,  Hall,  and  Rohsenow 
is  only  applicable  to  water,  where  the  new  method  developed  in  this 
study  Is  applicable  to  any  fluid. 

Latsh's  data  for  one-phase,  forced-convection  and  subcooled 
boiling  is  shown  In  Figure  9.  This  figure  shows  where  the 
calculated  incipient  boiling  point  is.  Notice  that  this  point  is 
where  the  data  makes  a  bend  upward.  This  upward  bend  is  also 
evident  in  the  computer  generated  lines. 

Wright's  data  for  low-quality  and  high-quality  flow  Is  shown  in 
Figures  10  through  13.  The  accuracy  of  the  low-quality  predictions 
clearly  shows  that  convective  evaporation  is  taking  place  under 
low-quality  conditions.  In  fact,  in  Appendix  A  where  a  sample 
calculation  is  shown,  forced-convective  evaporation  comprises  all 
of  the  heat  transfer  at  qualities  lower  than  five  percent  under  the 
flow  condition. 

It  should  be  noted  that  this  method  is  similar  to  that  develop¬ 
ed  by  Bjorge,  et.  al.  They  chose  X  =  0.05  as  the  boundary  between 
equation  (IV-23)  and  (IV-22)  and  did  not  use  equation  (V-l). 

However,  the  flow  visualation  experiment  done  in  this  study  and  a 

T231 

comparlsion  of  the  equations  with  Wright's1  J  data  showed  that  the 
boundary  between  equation  (IV-23)  and  (IV-22)  should  be  between 
zero  and  one  percent  quality.  For  computational  ease,  this  value 
was  chosen  as  zero  percent  quality.  Furthermore,  sample 
calculations  (see  Appendix  A)  showed  that  nucleate  boiling  is 


«UuU«' 


suppressed  at  relatively  low  qualities.  If  equation  (IV-22) 
continues  to  be  used  for  suppressed,  nucleate-boiling  conditions 
the  nucleate-boiling  term  is  actually  being  subtracted  off.  This, 
of  course,  is  inaccurate.  Therefore  equation  (V-l)  must  be  used 
under  suppressed,  nucleate-boiling  conditions. 


FLOW  STABILIZATION 

A  method  of  stabilizing  the  flow  is  to  prevent  slug  flow  from 
occurring.  This  can  be  don«  by  inducing  an  artif teal,  body  force 
so  that  stratified  flow  occurs  Instead  of  slug  flow.  The  simplest 
way  of  accomplishing  this  is  to  induce  swirl  in  the  flow  so  that 
there  is  a  large  centrifugal  force.  To  induce  swirl  a  twisted  tape 
inside  the  tube  can  be  used.  If  the  wicklng  height  of  the  fluid  is 
used  as  a  transition  point  between  slug  and  stratified  flow^: 


r  20 
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The  tangential  acceleration  due  to  fluid  swirl  is: 


(V-6) 


(V— 7) 
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boiling.  D»ta  from  Reference  21. 


2  bar 


bulk-fluid  temperature  versus  the  heat  flux 
subcooled  boiling.  Data  from  Reference  22. 


flux,  high-quality  boiling.  Data  from  Reference  23. 


where  u  is  the  tangential  velocity  at  the  tube  wall,  u  can  be 
tw  tw 

calculated  from: 


utv  *  fcT  r(RPM) 


(V-8) 


i 

: 


where : 

u  (V-9) 

rpm  «  15 

u  is  the  axial  velocity  and  Y  is  the  number  of  diameters  per 
a 

180-degree  twist  (twist  ratio).  Therefore: 


(V-10) 


Setting  equal  to  g: 

tang 


or 


0.31p«r 


The  Y  value  calculated  for  the  experiment  with  a  mass-flow  rate 
of  5  X10-3  kg/s  and  an  evaporator  diameter  of  0.25  inches  was  0.66. 
All  attempts  to  make  a  twisted  ribbon  with  this  twist  ratio,  with 
the  facilities  available,  failed.  Therefore,  this  method  could  not 
be  experimentally  verified. 


PRESSURE-DROP  CALCULATION 

r  24 1 

In  a  report  published  by  Bae  et  al.  the  pressure  gradient 
in  a  tube  for  two-phase  flow  is: 


dP  dp\  ,  dp\  dp\ 

*  HzJ,  .  +  +Tz ) 

/frict  /mom  /grav 


where  dP/dz)frict  is  the  frictional-pressure  gradient,  dP/dz)mojjj  is 

the  pressure  gradient  due  to  a  change  in  momentum  as  the  liquid 

changes  to  a  vapor,  and  dP/dz)  is  the  gravitational-pressure 

grav 

gradient.  The  frictional-pressure  gradient  can  be  calculated  with: 
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The  momentum-pressure  gradient  can  be  found  with: 


(V— 14) 


And  the  gravitational-pressure  gradient  can  be  found  with: 
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where  dy  is  the  differential  vertical  drop. 

These  equations  were  programmed  into  a  computer  routine  and 
compared  with  data  from  the  experimental  apparatus.  The  equations 
predicted  a  good,  average,  pressure  drop.  However,  the  accuracy  of 
the  equations  could  not  be  measured  because  of  the  severe  pressure 
fluctuations  encountered  in  the  test  section. 


VI  EVAPORATOR  DFSIGN  PROGRAM 


The  methods  outlined  in  section  V  were  used  to  develop  a 
program  to  design  straight-tube  evaporators  for  electronics  cooling 
at  10°C  using  either  ammonia  or  freon-11  or  at  100°C  using  water. 
Other  fluids  can  be  used,  but  require  the  addition  of  their  fluid 
properties  to  the  program.  The  user  of  the  program  chooses  a 
working  fluid,  and  supplies  the  total  power,  diameter  of  the  tube, 
inlet  temperature,  allowable  maximum  temperature  of  wall,  and  the 
exit  quality  or  the  mass-flow  rate.  The  program  then  calculates 
the  mass-flow  rate  needed  or  the  exit  quality,  whether  it  is 
possible  to  design  an  evaporator  with  the  given  inputs,  the  lengths 
of  the  one-phase,  subcooled,  and  two-phase  sections  and  the  total 
length,  the  wall -temperature  differences  for  each  section,  the 
average,  heat-transfer  coefficient  for  each  section,  the  pressure 
drop  for  each  section,  and  the  total  pressure  drop. 

To  calculate  the  length  of  an  evaporator,  the  following  methods 
can  be  used.  For  an  evaporator  subjected  to  a  constant  heat  flux, 
the  incipient-boiling  point  must  be  calculated  first  using  equation 
V-4.  Then  hop  and  htp  are  calculated  using  equations  IV-27,  IV-23, 
IV-24  and  V-l.  And  finally,  using  the  allowable  aT>  average  htp, 

Q,  and  d,  the  length  can  be  found.  For  an  evaporator  subjected  to 
a  constant  wall  temperature,  calculate  the  average  two-phase 
heat-transfer  coefficient  using  equations  IV-23,  IV-24,  and  V-l. 
Since  on-phase  heat-transfer  cannot  occur,  this  will  be  the  average 
heat-transfer  coefficient  for  the  entire  tube.  Next,  calculate  the 
heat  flux  using  the  given  wall  temperature.  And  finally,  calculate 
the  length.  The  pressure  drop  can  be  calculated  using  equation 


i 
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Example  1:  Water 

ENTER  FLUID:  WATER ( 1 > r AMMON IA ( 2 ) , F-l 1 ( 3  )  1 
ENTER  POWER  <W>  »1000 

ENTER  DIAMETER  (M)  »1.27E-2 

ENTER  INLET  TEMPERATURE  <K)  >>363.0 
ENTER  DELTA  T  <K>  >>10 

DO  YOU  WANT  TO  ENTER  EXIT  QUALITY  <1)  OR  FLOW  RATE  ( 

EXIT  QUALITY  FRACT ION  >>.8 

THE  MASS  FLOW  RATE  IS  0 . 5408E-03 ( ks/s ) 


INITIAL  FLOW  IS  LAMINAR 
BOILING  STARTS  IMMEDIATELY 


ONE-PHASE 

SUBCOOL  BOILING 

TWO-PHASE 

LENGTH  (m) 

0 . OOOOE+OO 

0 . 4014E-02 

0.1697 

DELTA  T  <  K ) 

0 . OOOOE+OO 

18.74 

10.00 

h  <W/m~2  K) 

233.6 

7703. 

0. 1443E+05 

TOTAL 

0.1737 


DELTA  P  <PSI) 


0.0000 E+00  0. 1441E-06 


0.3694E-02  0.3694E-0 


Example  2:  Ammonia 


ENTER  FLUID:  WATER < 1 ) , AMMON IA < 2 ) , F- 1 1 ( 3 ) 2 
ENTER  POWER  <W>  >>1000 

ENTER  DIAMETER  <M)  >>1.27E-2 

ENTER  INLET  TEMPERATURE  (K)  >>273 
ENTER  DELTA  T  <K)  >>10 

DO  YOU  WANT  TO  ENTER  EXIT  QUALITY  (1)  OR  FLOW  RATE  (2>?>>1 

EXIT  QUALITY  FRACTION  ».8 

THE  MASS  FLOW  RATE  IS  0 . 9731 E-03 ( kg/s ) 

INITIAL  FLOW  IS  LAMINAR 

BOILING  STARTS  IMMEDIATELY 

ONE-PHASE  SUBCOOL  BOILING  TWO-PHASE  TOTAL 

LENGTH  ( m )  O.OOOOE+OO  0.1626E-01  0.3433  0.3596 

DELTA  T  (K)  O.OOOOE+OO  10.93  10.00 

h  ( W/m"2  K)  100.3  6377.  6970. 


DELTA  P  <PSI> 


O.OOOOE+OO  0. 1352E-05 


0.1591E-02  0. 1593E-02 


Example  3:  Freon-11 


ENTER  FLUID:  WATER ( 1 > , AMMON IA< 2 ) , F-l 1 < 3 >3 
ENTER  POWER  (W)  >>1000 

ENTER  DIAMETER  <M)  »1.27E-2 

ENTER  INLET  TEMPERATURE  <K>  »273 
ENTER  DELTA  T  (K)  >>15 

DO  YOU  WANT  TO  ENTER  EXIT  QUALITY  (1)  OR  FLOW  RATE  <2>?»1 

EXIT  QUALITY  FRACTION  ».8 

THE  MASS  FLOW  RATE  IS  0 . 6377E-02 < kg/s > 

INITIAL  FLOW  IS  LAMINAR 

BOILING  STARTS  IMMEDIATELY 

ONE-PHASE  SUBCOOL  BOILING  TWO-PHASE  TOTAL 

LENGTH  <»>  0 . OOOOE+OO  0.1967E-01  0.3397  0.3593 

DELTA  T  (K)  0. OOOOE+OO  17.98  15.00 

h  ( W/m"2  K)  35.73  3880.  4650. 


DELTA  P  (PSI) 


0. OOOOE+OO  0 . 9424E-05 


0 . 8642E-01  0.8643E-01 


Exaaple  A:  Subcooled  Water 


ENTER  FLUID:  WATER < 1 ) , AMMONIA<  2) , F-l 1 ( 3 )  1 
ENTER  POWER  (W)  >>500 

ENTER  DIAMETER  (M)  »1.27E-2 

ENTER  INLET  TEMPERATURE  <K)  >>330 
ENTER  DELTA  T  <K>  >>15 

DO  YOU  WANT  TO  ENTER  EXIT  QUALITY  (1)  OR 
ENTER  MASS  FLOW  RATE  <kg/s)  »6.0E-3 

ENTIRE  SECTION  IS  SUBCOOLED 


VII  CONCLUSIONS 


1.  From  the  review  of  the  forced-convectinve  boiling  process: 

la.  The  gravitational  acceleration  has  very  little  effect  In 
the  subcool ed-bol ling  and  two-phase,  forced-covectlve  region. 
Therefore,  the  lack  of  a  gravitational  acceleration  (o-g)  should 
not  affect  the  calculation  of  local  heat-transfer  coefficients. 

lb.  The  gravitational  acceleration  does  affect  the  saturated 
nucleate-boiling  region.  However,  since  the  gravitational  effect 
tends  to  balance  out  around  the  circumference  of  the  tube  and  since 
this  region  is  generally  very  short,  the  use  of  equations  developed 
for  1-g  accelerations  can  be  used  to  calculate  the  local 
heat-transfer  coefficient  without  much  loss  of  accuracy. 

2.  From  the  flow-visualization  study,  placing  a  twisted-tape 
insert  into  the  evaporator  should  stabilize  the  flow  by  Inducing  a 
large  centrifugal  acceleration  which  prevents  slug  flow  from 
occuring.  The  twist  ratio  must  be  large  enough  to  ensure 
stratified  flow  under  all  conditions.  This  twist  ratio  can  be 
calculated  using  equation  V-ll. 

3.  The  length  of  a  straight-tube  evaporator  can  be  calculated  for 
a  constant  heat  flux  or  a  constant  wall  temperature  by  using  the 
methods  developed  in  this  study  and  summarized  in  Section  VI. 
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APPENDIX  A 


This  appendix  shows  how  the  calculations  were  done  to  obtain 
the  solid  lines  on  the  graphs.  More  specifically  sample 
calculations  are  shown  for  Figures  8,  9,  and  10.  In  addition,  a 
sample  twist  ratio  calculation  is  done. 


SAMPLE  CALCULATION:  Subcooled  Boiling,  Figure  8 


Data: 


Equation: 


Calculation: 


k  =  0.68  W/m  K 
p  =  958.0  kg/m3 

~  O 

pv  =  0.6  kg/m 

h .  =  2.258X106  J/kg 

fs  «;  2 
P  =  l.OlUXlO"7  N/m 


.I 

V  =  2.8X10  kg/m  s 
a  =  5.89X10-2  N/m 

d  =  1.27X10*2  m 
Cp  =  4. 21X10 3  J/kg  K 


,  0  23fo.66w/m K .Vail. 

\l.27X10’2m/\(2.2 

*(' 


27X10~2m) 


25  8X10W(  2 . 8X10" \g /ms ) 


4 . 21X10 3J/kgK)  ( 2 . 8X10~Stg/ms  )\°  ’ 69 


0.68¥/mK 


L) 


'( 1 . 0l4X105B/m2 )  ( 1. 27X10~~2- N  N  0’31 


^(1.014 


89X10"2N/m 


/958.0kg/m3  \°'31 

\  0.6kg/m3  “  V 


hNB  =  2.45q°’69  W/m2K 
4  *  hNBATsat 


q  =  18.07(30) 

£  P 
q  =  1.06X10°  W/m 


SAMPLE  CALCULATION:  Subcooled  Boiling ^  Figure  9 


One-Phase  Forced  Convection: 

Data:  T  -T„t  =60  K 

w  FL 

Tfl  =  30  C 
u  =  0.5  m/s 
d  »  1.2X10-2  m 


=  983.3  kg/m3 
y  =  4.71X10"14  kg/m  s 
k  =  0.654  W/m  K 
Pr  =  3.01 


Equation: 


hpc  =  0.23Ee°-V-1,i. 


Calculation: 


ref 


=  a  i2|32.  =  6o  c 


Re 


=  _(983.3kg/m3)(0.5m/s)(l.2X10~2m) 

_  u  ”  .  _  _ 


.  -i* 

4.  71X10  kg/m  s 


Re  =  1.253X10 


4 


hpc  =  0. 023(1. 253XloS0*8(3.0l)0,1* 


0.654w/mK 

1.2X10_2m 


hpc  =  3697  W/m2  K 

q  =  hpcAT 


q  =  (3697W/m2K)(60K)  =  2.22X105  W/m2 


Incipient  Boiling  Point: 

Data:  k^  =  0.68  W/m  K 

p„  *  958.0  kg/m3 

x,  s 

h.  =  2.258X10  j/kg 
Jg  _h 

y  =  2.8X10”  kg/m  s 


O  =  5.89X10"2  N/m 
Pr  =  1.66 

v  f,  =  1.6719  m3/kg 

*6  P 

d  =  1.2X10”^  m 


Equation: 


8oT 


AT 


sat, IB 


sat  ^ghFC 


k.h . 

I  fg 


(Pr)' 


Calculation: 


(958kg/m3)(0.$m/s)(l.2X10~2m) 
2.8X10"4kg/m  s 


Re  =  2.053X10 


hpc  =  0.023Re°-8Pr°*4  | 


h__  *  0 .023 (2.05 3X104 )° ' 8(  1.66 )° ’ 4  °'68W/m!J 

1.2X10  ra 


h__  =  1*498  W/m  K 

r  C 


AT 


sat , IB 


_  (8)(5.89Xl0~2N/m)(373K)(l.6719m3/kg) 
(0.68w/mK)(2.258xi06j/kg) 


X  (2*l*l*7W/m2K)(l.66)2 


AT  .  TTJ  ■  2.37  K 
sat, IB 


Subcooled  Boiling: 

Data:  p^.  =  0.6  kg/m3  P  =  2XlO^N/m2 

Cp  =  4.21X10 3  J/kg  K  AT  *  20  K 

sat 


68 


Equation: 


h_T_  = 
NB 


q(l.2X10~2m) 


=  0.25  (—  ---  - - - 

\  1. 2X10~  m  (2.25  8X10°  J/kg )  ( 2 . 8X10“\g/ms ) 


2 1X10  3 J/k  gK ) ( 2 . 8X10~^kg/ms )\° ' 69 


0.68W/mK 


X 


/ (2X10^N/m2)(l.2X10~2m) 
\  5.89X10_2N/m 


958kg/m3 

0.6kg/m3 


h 


NB 


3.02q 


0.69 


q  =  \bAT 


sat 


q  =  3.02q0-69AT 


sat 


q  =  35.4AT  3,23 

S  8-X 

q  =  35.4(20  K)3-23  =  5.65X105  W/m2 


SAMPLE  CALCULTION:  Forced-Convective 
Data:  k  =  0.68  W/m  K 

=  958  kg/m 3 
,  Pv  =  0.6  kg/m3 

D0  =  2.8X10-^  kg/m 


Boiling  Figure  10 

h  =  2.258X106  J/kg 
*  & 

Cp  =  4. 21X10 3  J/kg  K 
a  =  5.89X10“2N/m 
s  p  =  1.2TX10"5  kg/m  s 


a  =1.91X10  m 


q  =  1*3525 


P  =  1.U5X105  N/m2  Pr  -  1.66 
▼f  =  1.67  m3/kg  x  =  0.03 


Equations: 


0.69/n_.v0. 69/^0.31 


(¥■)  (?)  (?-■) 


hp^  *  0.023Re^‘ ^Pr^ * ^ ~  (one-phase) 


ReJl0‘9pr*F(Xtt)kJl 

qFC  * - T, ^ - ATsat  ( two~Phase ) 


rtv-°-»(Va-°t 


1  =  qFC  +  qjjg 


(■•(• 


sat. IB1 
AT  „  . 


-  V  ) 


Calculations: 


Re  „  md  =  ( ,  21kg/s  )  (1.91X10  m) 


Au  Wl.91X10"2m)2(2.8xl0"\s;g/ms) 


Re  =  5X10 


h  =  0. 023(5X10 ^)0-8(l. 66)0,1* 

1.91X10  m 


h__  =  5760  W/m2K 


-  (8)(5.89X10~2N/m)(373K)(l.67nrVkg) 
(0.68w/mK)(2.758X106J/kg) 


sat, IB 


X  (5760W/m2K)(l.66)2  *  2.71  K 


hNB  * 


U)(0.68w/aK) 
6  """ 


\0.69 

)/ 


=  0.25/^ 

\(2. 

x  / d.  21X103J/kgK)  (2. 8X!0~^kg/ms )\  0,69 

V  0.68W/mK  / 

v/(l.^5X10^N/m2)(l.91X10~2mA°  * 31 
\  5.89X10'2N/m  / 

X  iV 

\0 . 6kg/m3  / 

h  =  2.735<i0'69  *  2.735(^3525W/m2)0-69 


0.31 


hNB  =  U3Ul  W/m2lC 

X  -( Q-^kg/m3  \^*^/ 2.8X10JWms  N0*1  /l-0.03\°'9 
%t  \958kg/m 3)  \1.27X10_3kg/ms/  \  0.63  / 

Xtt  =  0.7788 

Re£  =  Re(l  -  x)  =  -U.85X105 

/  n  1  \0.32 

F(Xtt)  =  0.15  +  2.0  ypfwj 

F(Xtt)  =  0.5176 


F2  *  5(1.66)  +  51n(l  +  5(1.66))  +  2.51n(0.0031 


71 


since  AT  T_  =  2. 
sat,  IB 

Therefore: 

SAMPLE  CALCULATION: 
Data: 

Equation: 


X  ( U. 85x10 11 )° • 812 ) 


f2  *  26.91 


*FC 


_  ( U . 85X10^ )° ‘ ^ ( 1. 66 ) ( . 5176 ) ( . 068W/mK) 

/  A  1  ^  /  i  A  1  *•«*  \ 


(26.91) (1.91X10  m) 


%c  - 


9  -  k3525 


l87434Taat 


+  I^UIAT 


sat 


C 


2.7 

AT 


AT 


sat 


=  2.1+9  K 


71  >  ^Tsat,  toiling  is  suppresed  and  q  =  qF 
h  a  1871+3  W/m2K 

Twist  Ratio 

A  =  5X10"3  kg/s  d  =  0.25  in  *  6.35X10 

P£  *  955  kg/m3  a  =  5.89X10"2  N/m 


«  -  4  -  sno-y. _ 

a  pA  (955kg/m5)(u/l+(6. 35X10“  3m)2) 
u  -  1.65X10"1  m/s 


Calculation: 


Y  = 


(- 


. 3l)(955kg/Q3)(6. 35X10'  ' 


(2)(5.89X10‘2  H/m) 


WJ 


1.65X10_1m/s 


Y  =  0.66  diameters/l80  degree  twist 


APPENDIX  B 


This  appendix  is  a  listing  of  the  straight-tube  evaporator 
design  program.  It  is  explained  in  section  VI  and  several  sample 
runs  are  shown  on  pages  55-58. 


PROGRAM  EVAPD 

DIMENSION  X(25),HBH(25),DTBH(25),HTP<25) 

C 

-CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

C 

C  This  program  designs  an  evaporator  given  a  constant  heat  flux 
C 

CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 
C 

C  This  is  the  data  section 

DATA  GR, IMAX,FRC/9.8,20, .036/ 

CALL  ASS  IGN( 60 ,'TT:',3> 

CALL  ASS IGN<  6 ,  ' TT: ' ,3 ) 

WRITE  (60,7) 

7  FORMAT < 1H  , 'ENTER  FLUID:  WATER ( 1 > , AMMONIAC 2),  F -11(3)',$) 

READ  <6, A)  IFLD 
IF  < IFLD . EQ . 1 )  GO  TO  1 
IF  ( IFLD . EG . 2 )  GO  TO  2 
IF  <  IFLD . EQ . 3 )  GO  TO  3 
C 

1  CALL  H20(TS,PR, XMUL , XMUV , HFG , RHOL ,RHOU,CPL,SIG,VFGfXKL,P) 

GO  TO  6 

2  CALL  AM(TS,PR, XMUL , XMUV , HFG , RHOL , RHOV ,CPL,SIG,VFG,XKL,P) 

GO  TO  6 

3  CALL  F11(TS,PR, XMUL , XMUV , HFG,RHDL,RHQU,CPL,SIG,UFG,XKL,P) 

GO  TO  6 

C 

C  The  symbols  in  the  data  section  are:TS=saturation  temper ature , GR= 

C  gravitational  accel leration, PR=Prandtl  number ,XMUL~liquid  viscosity 
C  XMUV=vapor  viscosity , HFG^specif ic  heat  of  vapor ization,RHOL-liquid 
C  density ,RH0V=vapor  density ,CPL=1 iquid  specific  heat,SIG=surface 
C  tension, VFG-specific  volume, IMAX=number  of  two-phase  sections, XKL= 

C  liquid  conductivity , FFC=fr ict ion  coefficient 
C 

6  WRITE  (60,10) 


10 

FORMAT ( 1H  , 'ENTER  POWER 

READ  (6, A)  Q 

WRITE  (60,30) 

(W) 

V 

V 

30 

FORMAT ( 1H  , 'ENTER  DIAMETER 

READ  (6, A)  D 

WRITE  (60,60) 

(M) 

>>',$> 

60 

FORMAT ( 1H  , 'ENTER  INLET  TEMPERATURE  (K) 
READ  (6, A)  T I 

IF  (TI.GE.TS)  GO  TO  2120 

WRITE  (60,62) 

V 

V 

<e 

-62 

FORMAT ( 1H  , 'ENTER  DELTA  T 

READ  (6, A)  DELT 

WRITE  (60,65) 

(K) 

»'  ,$) 

65 

FORMAT ( 1H  , ' DO  YOU  WANT  TO  ENTER 

EXIT  QUALITY 

♦  FLOW  RATE  (2)?»/,t) 
READC6 , A )  ID 
IF  (ID.EQ.l)  GO  TO  01 


90  FORMAT ( 1H  , 'ENTER  MASS  FLOW  RATE  (kg/s)  »',*> 

READ (6,*)  XMD 

GO  TO  99 

91  WRITE  (60,95) 

95  FORMATdH  ,'EXIT  QUALITY  FRACTION  >>',♦) 

READ ( 6 , A  )  XE 

XMD=Q/(CPLA(TS-TI)+HFGAXE) 

WRITE  (60,96)  XMD 

96  FORMAT( 1H0, 'THE  MASS  FLOW  RATE  IS  ' , G10 . 4 , ' ( kg/s ) ' ) 

GO  TO  150 

99  QOP=XMDACPLA(TS-TI) 

IF  (QOP.GE.Q)  GO  TO  2100 
XE=(Q-QOP)/( HFGAXMD ) 

IF  (XE.GT.1.0)  GO  TO  2000 

WRITE  (60,101)  XE 

101  FORMAT< 1H0, 'THE  EXIT  QUALITY  IS  ',G10.4) 

C 

CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

C 

C  This  section  calculates  the  one- phase  Reynolds  number 

C 

C 

150  G=XMD/(3. 1416/4. 0A<DAA2) ) 

RE=GAD/XMUL 

1=2300 

C 

C  Calculate  the  one-phase  heat  transfer  coefficient 
C 

CALL  SHOP (RE,XKL,D,PR,HOP) 

C 

C  Find  the  incipient  bioling  piont 

C 

CALL  TIB(SIG,TS,VFG,XKL,HFG,HOP,DTIB) 

C 

C  Find  the  amount  of  energy  needed  to  reach  saturation 

C 

QSC=XMDACPLA(TS-TI) 

IF  (QSC.GE.Q)  GO  TO  2020 
C 

C  This  part  divides  the  two-phase  section  into  parts 
C 

QTP=Q-QSC 
XMAX=1 .0 

X( IMAX ) =QTP/ (XMDAHFG ) 

IF(X( IMAX) .GE.XMAX)  GO  TO  2000 
IMPO* IMAX+1 
DELX=X ( IMAX ) /IMAX 
X(1)=0.0 
DO  300  1=2, IMPO 
X(I)«X(  I-D+DELX 
300  CONTINUE 

C 

C  Calculate  two-phase  heat  transfer  coefficients 
C 


C  Estimate  the  length 
C 

CALL  SHMN< XKL,D,Q,XMUL,HFG, DELI, PR, RHOV, RHOL, P,SIG,HMN) 
ZLTG*Q/(3. 1416ADAHMNADELT > 

— -350  OOA*Q/<3. 1416ADAZLTG) 

C 

C  Begin  two-phase  calculations 
C 

400  HTPA*0. 0 

DO  500  1*1 , IMPO 

CALL  SHBH(DTIB,TS,X< I),XIA,XKL,D,Q,XMUL,XMUV,HFG,ZLTG, 

A  PR , RHOL , RHOV ,P,SIG,HOP,RE,HBH(  I)  > 

HTPA*HIPA+HBH< I) 

500  CONTINUE 

HTPA=HTPA/ IMPO 
C 

C  Calculate  lengths 
C 

ZLT=Q/(3. 1416ADAHTPAADELT ) 

Q0A=Q/<3. 1 416ADAZLT ) 

IE  (ABS(ZLT-ZLTG) .LE.0.01 >  GO  TO  510 

ZLTG=ZLT 

GO  TO  400 

510  DTOP=QOA/HOP 

ZLOP=CPLAXMDA<  DT IB+TS-QOA/HOP-T I ) / ( Q0AA3 . 1416AD ) 

IF  <  ZLOP . LE . 0 . 0 )  GO  TO  2040 

ZLSCB*  <  TS- <  DT IB +TS-QO A/HOP ) ACPLAXMD ) / <  Q0AA3 . 1 41  GAD) 
ZTEST*ZLOP+ZLSCB 
IF  (ZTEST.GE.ZLT)  GO  TO  535 
GO  TO  530 

520  ZLSCB=(TS-TI)AXMDACPL/<  Q0AA3. 1 41 6AD ) 

ZTEST=ZLOP+ZLSCB 
IF  (ZTEST.GE.ZLT)  GO  TO  535 
530  ZL TP* ZLT-ZLOP- ZLSCB 

GO  TO  540 
535  ZLTP=0 . 0 

ZLSCB=ZLT-ZLOP 


C 

C  Calculate  delta  wall  temperature 
C 

540  DTSCB«Q0A/HBH<1 ) 

DTTP*QOA/HTPA 

C 

C  Calculate  pressure  drops 
C 

CALL  SDPOP (XMD, RHOL, D,XMUL, ZLOP ,FRC, RE, T , DPOP , DPOPE > 

CALL  SDPOP ( XMD , RHOL ,  D , XMUL , ZLSCB, FRC, RE, T,DPSC,DPSCE) 
DELL*ZLTP/IMAX 

CALL  SDPTP(XMUL,XMUV,X, RHOL, RHOV, XMD, D, IMPO, VFG, DELL, DPIP) 
DPTE-DPOPE+DPSCE+DPTP 
C 

CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA^ 

-C 

C  This  section  prints  out  the  calculated  values 
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CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA* 

C 

WRITE  <60,1 000 ) 

1000  FORMATdHO, T20, 'ONE-PHASE' ,T32, 'SUBCOOL  BO IL ING ' , T50 , 

+  'TWO-PHASE' ,T62, 'TOTAL' ) 

WRITE  (60,1020)  ZLOP , ZLSCB , ZLTP , ZLT 

1020  FORMAT <1H0, 'LENGTH  < » > ' , T20 , G10 . 4 , T32 , G10 . 4 , T50 , G10 . 4 , T62 , G 10 . 4 ) 

WRITE  (60,1040)  BTOP, DTSCB , DTTP 

1040  FORMATdHO, 'DELTA  T  <  K )  '  ,  T20  ,  G 10 . 4  ,  T32 ,  G10 . 4  ,  T50 ,  G 10 . 4  ) 

WRITE  (60,1060)  HOP , HBH < 1 ) , HTPA 

1060  FORMATdHO, 'h  ( W/m''2  K ) '  ,  T20 , G1 0 . 4 ,T32 ,  G10 . 4 , ISO, G10 . 4  ) 

WRITE  (60,1080)  DPOPE , DPSCE , DPTP , DPTE 
1080  FORMATdHO, 'DELTA  P  (  PS  I  >  '  ,  T20 ,  G10 . 4  ,  T32  ,  G10 . 4  , 

+  T50,G10.4,T62,G10.4) 

GO  TO  3000 

2000  WRITE  (60,2015) 

2015  FORMATdHO, 'BURN  OUT  OF  EVAPERATOR') 

GO  TO  3000 

2020  WRITE  (60,2030) 

2030  FORMATdHO, 'THE  ENTIRE  SECTION  IS  ONE-PHASE  FLOW') 

GO  TO  3000 

2040  WRITE  (60,2050) 

2050  FORMATdHO, 'BOILING  STARTS  IMMEDIATELY') 

ZL0P=0 . 0 
DT0P=0.0 
GO  TO  520 

2100  WRITE  (60,2110) 

2110  FORMAT  (1H0, 'ENTIRE  SECTION  IS  SUBCOOLED') 

GO  TO  3000 


2120  WRITE  (60,2130) 

2130  FORMAT  (1H0, 'ENTIRE  SECTION  IS  SUPERHEATED  VAPOR') 

3000  STOP 

END 
C 

C A A AAAAAAAAAAAAAAAAAAA AAA A AAA A A A A AAAAAA A A AAAAAAAAAAAAAAAA A A A A A A A A A 
C 

C  Start  of  subroutines 
C 

CAAAAAAAAAAAAAA AAAAAAAAAAAAAAAA AAAAAAAAAAAAAAAAAAA AAAAAA AAA A AAAAAA 

-C  - 

SUBROUTINE  T IB ( S IG , TS , VFG , XKL , HFG , HOP , DT IB ) 

C 

C  This  calaculates  the  incipient  boiling  piont 
C 

DT  IB*8. OASIGATSAVFGAHOP/ <  XKLAHFG ) APRAA2 

-  -RETURN 

END 

C 

SUBROUTINE  SHOP(RE,XKL,D,PR,HOP) 

C 

C  This  subroutine  calculates  the  one-phase  heat  transfer  coefficient 
C 


T=2300.0 

IF  (RE.LE.T)GO  TO  120 
C 

C  This  section  calculates  the  turbulent  one-phase  heat  transfer 
C  coefficient 

--C 

WRITE  <60,100) 

100  FORMAT (1H0,' INITIAL  FLOW  IS  TURBULENT ' ) 

H0P=0.023AXKL/DA(REAA0.8)A<PRAA0.4) 

GO  TO  150 
C 

C  This  section  calculates  the  laminar  one-phase  heat  transfer 
C  coefficient 
C 

120  WRITE  (60,130) 

130  FORMAT < 1 HO, 23H INITIAL  FLOW  IS  LAMINAR) 

H0P=48.0AXKL/< 11 .OAD) 

150  RETURN 

END 


e 


c 

SUBROUTINE  SDPOP < XMD , RHOL , D , XMUL , ZLOP , FRC , RE , T , DPOP , DPOPE ) 

C 

C  This  section  calculates  the  one-phase  pressure  drop 
C 

V0P-XMD/<RH0LA3. 1416/4.0 AD AA 2) 

IF  (RE.GT.T)GO  TO  180 
DP0P=32 . 0AXMULAZL0PAV0P/DAA2 
GO  TO  190 

180  DPQP=FRCAZL0P/DAV0PAA2ARH0L/2 . 0 

190  DPOPE=DPOPA1 .45E-4 

RETURN 
END 
C 

SUBROUTINE  SHMN < XKL , D, Q , XMUL , HFG , DELT, PR , RHOU , RHOL , P, S IG , HMN ) 
C 

C  This  caculates  the  two-phase  heat  transfer  coefficient  using 
C  McNelly's  method 
C 

HMN=0.225AXKL/DA<  < DELTAD/ < XMULAHFG >  >AA0.69> 
HMN=HMNAPRAA0.69A<RH0L/RH0M-1 .0>AA0.31 
HMN=(HMNA(PA6.895E3AD/SIG) AA0.31 )AA3.226 
RETURN 
END 
C 

SUBROUTINE  SHBH< DT IB , TS , X , X IA , XKL , D , G , XMUL , XMUV, HFG, ZLTG , PR, 

+  RHOL  f RHOV , P , S IG , HOP , RE , HBH ) 

C 

-G^his  section  calculates  the  equation  derived  by  Bjorge,  Hall,  and 
C  Rohsenow 
C 

C  Initial  calculations 
C 

ZER0=0. 0 
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HB=0 . 225AXKL/0A < (0/(3. 141 6AXMULAHFGAZLTG > ) A AO. 69) 
HB=HBAPRAA0.69A<RH0L/RH00-1.0)AA0.31 
HB=HBA ( PA6. 895E3AD/S IG)AA0.31 
Q0A=Q/<3.1416ADAZLTG) 

AC=0 .01 
C 

X INC=10000 . 0 
DTG=X INC 
C 

C  Test  for  flow  regime 
C 

IF  (X.GT.0.0)  GO  TO  450 
C 

C  Subcooled  Regime 
C 

410  E=(1.0-(DTIB/DTG)AA3)AA2AHBAA2 

C=H0PAA2ADTGAA2 
A=Q0AAA2 . O-C 
C 

C  Test  for  negative  square  root 
C 

IF  ( A . LE . ZERO )  GO  TO  430 
DTBH=AAA0.5/(EAA0.5) 

ER=DTBH-DTG 

IF  <ER.LT.ZER0>  GO  TO  430 
ERAB=ABS ( ER ) 

C 

C  Test  for  accuracy 
C 

IF  <ERAB.LE.AO  GO  TO  440 
DTG=DTG+X INC 
GO  TO  410 

430  XINOXINC/IO.O 

DTG=DTG-XINCA10.0+XINC 
GO  TO  410 

440  HBH=QOA/DTBH 

GO  TO  490 
C 

C  High  quality  section 
C 

450  REL=REA  < 1 . 0-X ) 

F0=50.0 
FT  =  1 1 25 . 0 
C 

C  Test  for  proper  equation  use 
C 

IF  (REL.GE.FT)  GO  TO  470 

IF  (REL.GE.FO)  GO  TO  460 

FTW0*0. 0707APRARELAA0 . 5 
GO  TO  475 

460  C*1.0+PRA<0.0964ARELAA0. 585-1 .0) 

FTW0=5.0APR+5.0AAL0G(C) 

GO  TO  475 


470 


FIW0=5.0APR+5.0AAL0G<1 .0+5.0APR) 
FTW0=FTW0+2.5AAL0G<0.0031ARELAA0.812) 

475  XTT=  <  RHOL/RHOV )AA0.5A( XMUV/XMUL >  AAO .  1 

XIT=XTTA(X/< 1 .0-X) )AA0.9 
FXTT*0.15A<XTI+2.0AXTTAA0.32> 
HFC=RELAA0.9APRAFXTTAXKL/<FTWDAD> 

480  DTBH= (QOA-HBADTGA ( 1.0-<DTIB/DTG)AA3))/HFC 

ER=DTBH-DTG 

IF  (ER. LI. ZERO)  GO  10  485 
C 

C  test  for  accuracy 
C 

ERAB=ABS(ER> 

IF  (ERAS. LE. AC)  GO  TO  487 
C 

C  Test  for  negative  difference 
C 

IF  (ER.LE.ZERO)  GO  TO  485 

DTG=DTG+XINC 

GO  TO  480 

485  XINC=XINC/10.0 

DT6=DTG-X  INCA1 0 . 0+X INC 
GO  TO  480 

487  IF  (DTBH.LT.DTIB)  GO  TO  488 
HBH=QOA/DTBH 

GO  TO  490 

488  HBH=HFC 
DTBH=QOA/HBH 

490  RETURN 

END 
C 

SUBROUTINE  SDPTP ( XMUL , XMUV, X , RHOL , RHOV T XMD , D , IMPO , VFG , DELL , 
+  DPTP ) 

DIMENSION  X ( 25 ) 

C 

CThis  section  calculates  the  two-phase  pressure  drop 
C 

DPF=0 . 0 

DO  920  1=2, IMPO 

A=5.7A(XMUL/XMU0)AA0.0523A(1.0-X< I) )AA0.47 
A=AAX (  I)AA1.33A(RH0V/RH0L)AA0.261 
B=8.11A<XMUL/XMUV)AA0.105A<1.0-X< I) >AA0.94 
B=BAX (  I) AA0.86A<RHOV/RHOL)AA0.522 
C=X ( I ) AA1 . 8+A+B 

DP=0.1459AXMDAA2/(RH0VADAA5)ACA1 .347E-5 
DP=DPA< 1 .2732AXMD/(DAXMUV) > A A ( -0 . 2 ) ADELLA3 . 28 1 
DPF=DPF+BP 
*20 - -CONTINUE 

DPM=1.6211AXMDAA2/DAA4AX< IMPO ) AVFGA1 . 45E-4 

DPTP=DPF+DPM 

RETURN 

END 


SUBROUTINE  H20 ( TS , PR , XMUL f XMUV , HF6 , RHOL , RHOV , CPL , S  IG , VFG , XKL , P > 

TS=373 .15 

PR=1 . 66 

XMUL=0.28E-3 

XMUV= 1 . 27E-5 

HFG=2 . 258E6 

RH0L=958 . 0 

RH0V=0 . 6 

CPL=421 0 . 0 

S  IG=5 . 89E-2 

VFG=1 .6719 

XKL=0 .68 

P=15 . 7 

RETURN 

END 

SUBROUTINE  AM ( TS , PR , XMUL , XMUV , HFG , RHOL , RHOV, CPL , S IG f VFG , XKL f P > 

TS=283 . 0 

PR=3 .74 

XMUL=0.235E~3 

XMUV=0.965E-5 

HFG=1 . 2265E6 

RH0L=624 . 8 

RH0V=4 . 864 

CPL=4647.3 

S IG=2 . 307E-2 

VFG=0 . 206 

XKL=0 . 292 

P=89 .19 

RETURN 

END 

SUBROUTINE  FI  1 ( TS , PR , XMUL , XMUV , HFG , RHOL , RHOU , CPL , S IG , UFG , XKL „ P ) 

TS=283 . 0 

PR=4. 126 

XMUL=0 . 5E-3 

XMUV= 1 . 04E-5 

HFG=1 . 853E5 

RH0LS151 1.0 

RH0V=3 . 632 

CPL=858 . 3 

S IG=2 . 05E-2 

VFG= . 2747 

XKL*0.104 

P-8 .783 

RETURN 
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